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Abstract: We consider lepton pair production at a hadron collider in a class of 
effective theories with the relevant operators being four-fermion contact interaction. 
Despite the nonrenormalizable nature of the interaction, we explicitly demonstrate 
that calculating QCD corrections is both possible and meaningful. Calculating the 
corrections for various differential distributions, we show that these can be substantial 
and significantly different from those within the SM. Furthermore, the corrections 
have a very distinctive flavour dependence. And finally, the scale dependence of the 
cross sections are greatly reduced once the NLO corrections are taken into account. 
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1. Introduction 

While the Standard Model (SM) remains a very consistent explanation for nearly 
all data pertaining to high energy physics experiments, a few small discrepancies 
persist. Furthermore, there are theoretical issues that cannot even be addressed 
within the framework of the SM alone. Examples include the replication of the 
fermion families, the naturalness problem associated with the Higgs scale, charge 
quantization, the baryon asymmetry in the universe, the presence of dark matter 
etc.. Clearly, an answer to such vital questions may be obtained only in a model 
much more ambitious than the SM. Candidates for the role include, amongst others, 
supersymmetry [1], grand unification [2,3] (with or without supersymmetry) , family 
symmetries (gauged or otherwise) and compositeness for quarks and leptons [4]. In 
general, each such scenario (with its peculiar strengths and weaknesses) is associated 
with an individual set of tell-tale signatures. On the other hand, if the SM is to be 
a valid effective low-energy description of such bigger structures, one should be able 
to construct, within the ambit of the SM, operators that would encapsulate a class 
of remnant effects that could pertain to any of these scenarios. We illustrate this 
explicitly in the context of one of the above-mentioned scenarios. 

The replication of fermion families suggests the possibility of quark-lepton com- 
positeness. In such theories, the fundamental constituents, very often termed preons 
[5], experience an hitherto unknown force on account of an asymptotically free but 
confining gauge interaction [6]. At a characteristic scale A, this interaction would 
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become very strong leading to bound states (composites) which are to be identified 
as quarks and leptons. In most such models [7,8], quarks and leptons share at least 
some common constituents. Since the confining force mediates interactions between 
such constituents, it stands to reason that these, in turn, would lead to interactions 
between quarks and leptons that go beyond those existing within the SM. Well below 
the scale A, such interactions would likely be manifested through an effective four 
fermion contact interaction [9] term that is an invariant under the SM gauge group. 
A convenient and general parametrization of such interactions is given by [4, 10] 



r - — 

A 2 



% (q 7 M Pi q) (h, Pj I) + & (qPi q) (£Pj 



1.1) 



where i,j = L,R and Pi are the chirality projection operators. Note that the La- 
grangian of eqn.( |l.l| ) is by no means a comprehensive one and similar operators 
involving the quarks alone (or the leptons alone) would also exist. However, for our 
purpose, it would suffice to consider only eqn. flPI ). Within this limited sphere of 
applicability, the strength of the interaction may be entirely absorbed in the scale A, 
and the couplings rjij and £y canonically normalized to ±1. 

While we have sought to motivate eqn. (|l.l| ) in the context of compositeness, 
these are by no means the only scenarios ones that can give rise to such an effec- 
tive interaction lagrangian. As is well known, a four-fermion process mediated by a 
particle with a mass significantly higher than the energy transfer can be well approxi- 
mated by a contact interaction [9] term with a generic form as in eqn.( |l.l| ). Examples 
include theories with extended gauge sectors, leptoquarks [11], sfermion exchange in 
a supersymmetric theory with broken imparity [12] etc.. In all such cases, on inte- 
grating out fields with masses M, ^ A [13], a series of such higher- dimensional terms 
obtain. Those in eqn. (|l.l|) are just the lowest order (in A" 1 ) ones. 

Several points are in order here 

• In general, integrating out the heavy fields would result in an almost infinite 
number of higher-dimensional operators. The terms in eqn.( |l.l| ) are just some 
of the lowest order (in A" 1 ) ones relevant to four-fermion processes. 

• For a given model, the couplings rjij and £y generated by the process of in- 
tegrating out heavy fields would be related to each other. Such relations are 
model-specific and determined, to a large extent, by the flavour structure of 
the parent theory. As already indicated, we shall not consider any such flavour 
structure, but hold rj,£= ±1. 

• Even a requirement such as 577(2) ®[/(l) invariance for the effective Lagrangian 
would imply a relation between such terms, but involving different fields. How- 
ever, we shall concern ourselves with only those involving qq£ + £~ , noting that 
the results would be essentially the same for its 577(2) cousins. 
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• Low energy observables (for example, meson decays) lead to severe constraints 
[14] on several of these couplings and even more so, on their products. Al- 
though many of these bounds were derived in the context of specific ultraviolet 
completions, it is easy to see that they are equally applicable to the generic 
contact interactions. 

• Note that the vector-axial vector (VA) or scalar-psedoscalar (SP) nature of 
the 77- and ^-couplings do not necessarily reflect the spin of the integrated out 
field that led to such terms. 

Clearly, operators such as these could, in principle, lead to significant phe- 
nomenological consequences in collider experiments, whether e + e~ [18], eP [17] or 
hadronic. Given the higher-dimensional nature of C, it is obvious that the conse- 
quent effects would be more pronounced at higher energies. In other words, the 
fractional deviation over the SM expectations would be concentrated more at higher 
invariant masses M [15], with possibly some nontrivial dependence on the rapidity 
y as well [16]. For example, composite quarks and electrons have been proposed as 
a possible explanation for the high-Q 2 anomaly at HERA [17]. Some of the best 
constraints on compositeness, for example, came from the OPAL [18] and CDF [19] 
experiments. More recently the measurement of the Drell-Yan cross section [20] at 
high invariant masses set the most stringent limits on contact interactions of the type 
given in eqn.( |l.l|) . For example, within the V^A-type interaction scenario, the scale 
A is constrained to be A ^ 3.3-6.1 TeV [21,22], with the bound depending on the 
chirality structure of the operator. 

As is well known, QCD corrections can alter quite significantly the cross sections 
at a hadronic collider. Thus, these may have serious bearing on the discovery poten- 
tial of such experiments. Even for as simple a process as Drell-Yan, the leading order 
(LO) results seriously underestimate the cross sections. This has led to the incorpo- 
ration of the next-to- leading order (NLO) or next-to- leading log (NLL) [24,25] results 
in Monte Carlos codes [24] or event generators such as JETRAD [23]. However, no 
calculations exist for the higher order QCD corrections to cross sections mediated 
by a generic contact interaction. Consequently, all extant collider studies of contact 
interaction have either been based on just the tree level calculations, or, in some 
cases, have imlicitly assumed that the higher order corrections are exactly the same 
as in the SM. Clearly, this is an unsatisfactory state of affairs and, in this paper, we 
aim to rectify this by calculating the next-to-leading order QCD corrections for both 
the VA-ty\)e and the S'P-type contact interactions. 

It might be argued that, such theories being nonrenormalizable, any higher- 
loop calculation is fraught with danger. However, the very structure of such terms 
(namely the current-current form of the Lagrangian) along with the fact that only 
one of the currents comprises coloured fields allows us to reliably calculate QCD 
corrections. This holds not only for the specific interaction in question, but also for 
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other theories that satisfy the abovementioned criterion [26]. On the other hand, 
were we to attempt to calculate the NLO electroweak corrections, it is by no means 
certain that similar levels of reliability or usefulness could be reached. 

The rest of the article is organised as follows. In Section 2, we start by outlining 
the general methodology and follow it up with the explicit calculation of the NLO 
corrections to the differential distribution in the dilepton invariant mass. In the 
following section, we consider the rapidity distributions. Section 4 contains our 
numerical results. And finally, we summarize in Section 5. 



2. NLO corrections 

We consider lepton pair production at a hadron collider in the context of a generic 
contact interaction as exemplified by eqn.(|l.l|). In other words, the process is 

P( Pl ) + ( P\p 2 ) - l+(h) + l~{l 2 ) + X( Px ) (2.1) 

where Pi denote the momenta of the incoming hadrons and Zj those for the outgoing 
leptons. Similarly, the inclusive hadronic state denoted by X carries momentum 
px- The hadronic cross section is defined in terms of the partonic cross sections 
convoluted with the appropriate parton distribution functions /<f (#) and is given by 

da PlP2 f 1 f 1 f 1 da ab 

2S 7ra = dx i dx 2 f Pl {xx) f P2 (x 2 ) / dz 2 s —— 8(t - zx x x 2 ) 

«v , „x„Jo Jo Jo "V 

ab=q,q,g 

(2.2) 

with Xi being the fraction of the initial state hadron's momentum carried by the 
parton in question. In other words, the parton momenta fcj are given by ki = XiPi. 
The other variables are defined as 



S EE (pi + p 2 f S=(h + k 2 f Q 2 EE (h + l 2 f 

Q 2 Q 2 (2 ' 3) 

T EE — Z EE — — T = ZX\X 2 . 

S s 

It is convenient to symbolically cast the matrix element for the process as a sum of 
several current-current pieces with a "propagator" in between. In other words, 

M Total = ^ Had ■ p j ■ ^ Lept ( 2 - 4 ) 

3 

where the dots (•) denote Lorentz index contractions as appropriate and the propa- 
gators Pj are 

P-, = g- 2 !3^ = <J„A Pz = Q2_ M 2^ iMz T z " 9 "" Pz 

p - 4I ~p p - 4 * 

~ K 2 = K 2 = 
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With this definition, the partonic cross section for the process a{k\) +6(^2) —* j{q) + 



= ^ E / dP ^ i>f o6 ^'r ■ p s m ■ p*>(q 2 ) ■ v^ 11 ' , (2.6) 



where |.M afe ~^"'| 2 denotes the square of the hadronic current, and dPS m+ i the (m + 
l)-body phase space element, viz. 

dPS m+1 = JUT (<& 27r5+(pf)) 27r^(g 2 - Q 2 ) ^ ^ 

where n is the dimension of spacetime and S + (x) carries its usual meaning. The 
leptonic tensor, given by 

= |jjQ^2^ ( ^ (2.8) 
is straightforward to compute and leads to 

C^^{(-^ + W) C ' AQ2) (2.9) 

CspW) U=f = SP) 



with 



£^(Q 2 ) = yQ 2 , Czz{Q 2 ) = f ((V) 2 + (V)' )Q 2 

£ 7 z(Q 2 ) = ^^Q 2 , £ 5P (Q 2 ) = Q 2 



(2.10) 



In the above, a denotes the fine structure constant, while g%' R parametrize the cou- 
plings of the left- and right-chiral fermionic fields to the Z, viz. 

9a = 2 is a +9a) > fl'a = - e a tan 6W , flff = -2 T a 3 csc 29 w - e a tan 9 W 

(2.11) 

in terms of the Weinberg angle (6w) an d the electric charge (e a ) of the fermion in 
question. 

On substituting for Cjj/^u> in eqn.(|2.6|), we have, for the hadronic cross section, 



2^ r (r,Q 2 ) = - PAQ 2 )P;>{Q 2 )£AQ 2 )W^{t,Q 2 ) (2.12) 



jj'=7,.Z,VA,SP 
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where the hadronic structure function W is defined to be 




where Tjj depends upon the spin of the current in question, viz 



Tsp(q) = 1 



-9^ + 



Q 2 



) (j,j' = l ,ZZ ,VA) 



(2.15) 



To compute the Q 2 distribution of the dilepton pair, one then has to calculate the 
square of the hadronic matrix element \Ai ab ^^'\ 2 Tjji(q), preferably in a suitable 
frame so as to render the integrations over the phase space and z easy. 

Note that, the bare partonic coefficient function A is a singular object, suffering 
from each of ultraviolet, soft and collinear divergences. To handle these, we adopt 
dimensional regularisation. The renormalization procedure for the VA-type interac- 
tions is quite established and may be found, for example, in Ref. [26]. Note that, for 
the SP type interaction, one-loop corrections results in an extra term — proportional 
to ln^Q 2 //i 2 j — as compared to the VA interactions [29]. This, of course, is not 
unexpected, as contrary to the usual conserved vector currents, a scalar current is 
renormalized by QCD interactions. It is easy to see that this extra term is pre- 
cisely the one that is absorbed into the bare contact interaction coupling constant in 
defining the renormalized coupling 

To the ultraviolet regularized (and renormalized) operator, we must add the 
contribution from the real emission diagrams (gluon bremsstrahlung as well as the 
Compton process), and this exercise leaves us with a quantity that suffers only from 
collinear singularities. The latter, of course, can be removed through mass factoriza- 
tion. If fip be the factorization scale, then Drell-Yan coefficient functions, after mass 
factorization by A J J b , are related to the bare functions through 



A^(z,Q 2 ,e) = J2 r ca(z, » 2 F ,l/e) ®V db {z^ 2 F , 1/e) ® A%(z, Q 2 , /4) (2.16) 



c.d 



with the convolution defined to be 




(2.17) 



6 



The kernels T a b are related to the leading order Altarelli-Parisi splitting functions [27] 



Pfjj/(z) through 



r o6 (*,/4, l/e) = S ab 5(1 -z) + r§{z,fi) 

,2 \ e/2 



a s (/i 2 ) 



4tt 

Expanding eqn.( |2.l6D to order a s we have 



— ^qq ^ a s ® ^ a sA q ^ ^ 



(2.19) 



thereby leading us to the finite coefficient functions. The physical hadronic cross 
section may be obtained by folding these finite coefficient functions with appro- 
priate parton distribution functions. For the sake of completeness, we present the 
results below. To begin with, we denote the renormalized parton-parton fluxes by 
H ab (xi,x 2 ,/i F ) where 

H gq ( Xl ,x 2 ,» F ) = f^(x 1 , t 4)(f[ 2 (x2,^) + f q P2 (x 2 , t 4)) (2- 2 °) 

H qg (x 1 ,X 2 ,fi F ) = Hgq^Xx^Hp) . 

Then, the inclusive differential cross section may be expressed as 

da PlP2 v-^ f 1 f 1 f 1 

25— — — — (r, Q 2 ) = 22 / / dx 2 / dz 5(t - zx\x 2 ) 

dQ Jo Jo Jo 

[FsM+VA,q QsM+VA,q + ^SP,q QsP,q] 
GsM+VA,q = H g q(x 1 ,X 2 ,fi 2 F )^A q °q ) ' SM (z,Q 2 ,fI 2 F ) + d s A^'^ (z, Q 2 , J 

+ ^H gg (x 1 ,x 2 ,n 2 F ) + H gq (x 1 ,x 2 ,/j, F )ja s A$' SM (z,ii 2 F ) 

Qsp, q = H q g(x 1 ,x 2 ,fi 2 F ){A q f SP (z,Q 2 ,fi 2 F ) + a s A ( q f SP (z,Q 2 ,fi 2 F )} 
+ ^H qg (x 1 ,x 2 ,fx F ) + H gq (x 1 ,x 2 , ^ 2 F )^a s A { qg ) ' SP (z, fi F ) 

(2.21) 
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with the constants J-"sM+vA,q and T^ p containing all the dependences on the coupling 



constants and propagators, namely. 

,2 

q5 ~2e q g[ y q ~ v ^ 2 



? SM+VA,q 



v ^ Vr7 Q 2 -M. 



+\ (V ) 2 + (ah 2 ) ((g?f + (ai) 2 ) z Q 

{-e, J2 Vij + Zq(Q 2 - Mf) J2 Vijgl 

i,j=L,R i,j=L,R 



aA 2 



a 2 A 4 



i,j=L,R 



(2.22) 



i,j=L,R 



Q 2 

Zq= {Q 2 -M 2 ) 2 + T 2 z M 2 

For the vector-axial vector couplings, the results for the coefficient functions are 
analogous to the case of the SM [26] , namely 



AST =£*<i-*) 



N 



(2.23) 



11 



^[{-4 + 2 C (2)} ff (l-z)-(l + z)ln^-2 
{ (dfc + § *(1 " *) " (1 + *)} In ($) + 4(^) + " 



ln(z) 
1 - z 



_ 2tt 



2|l-2^ + 22 2 | In 



Zjl 2 F 



+ l + 6z-7z' 



(2.24) 



For the scalar-pseudoscalar couplings, on the other hand, the LO coefficient function 



is given by 



AST = 



n N ~\ > 

while at the next-to-leading order coefficient functions are 



(2.25) 



a(1),SP ^_ 
99 



V (1),SP 



4ttCV 
N 



2nT f 
N 



/ln(l-z) 



{-2 + 4C(2)}<5(l-z)+2(l-z) + 4(l + z 2 )(^ 



+2 



1 + z' 

JT^z) 



-In 



2 (l -2^ + 2^ 2 ) In 



' Q 2 (l-z) 



+ (1 -z)(7z-3) 



(2.26) 
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The SU(N) color factors in the above equations are 

N 2 -1 1 
Cf = ^iT> ° A = N > Tf= 2- (2 ' 27) 

3. Differential cross sections with respect to dilepton rapidity 

Having considered, in the previous section, the differential distributions with respect 
to the dilepton invariant mass, we now consider a second variable of interest, namely 
the rapidity of the pair. The latter can be expressed as 

Y = - log (^)= ~ log ( ^MILZ«cmh\ (31) 

with the second equality valid in the center of mass frame of the hadrons. Thus, the 
rapidity distribution may be computed simply by introducing the identity 



J d Ys{Y-\l 0t (^)) = l, 



Pi ■ q 



in eqn.(p.2|). This leads to 



2S ^dY (T ' F ' Q2) = 27 £ P i' m C "' m ^f- (T ' F ' Q2) - 

V j,j'=-y,Z,VA,SP 

(3.2) 

where the hadronic structure functions are given by 



dY 



t,Y,Q 2 ) = / dxx I dx 2 fa^xi) f P2 (x 2 ) / dzd(r - zx x x 2 ) 

„ u „■ „■/ >>0 JO JO 

J dPS m+1 \M ab ^"'\ 2 T jf (q)s(Y - ~ log ( 



P2 ■ q 



.pi ■ q 

(3.3) 

We start with the leading order case which involves just the calculation of the square 
of the matrix element for the process a{k\) + b(k 2 ) — > j{q)- The relevant phase space 
element corresponds to that for a (0 + l)-body final state, and 

J dPS 0+1 J dzS^Y -^log (^—^yjS(r - zxix 2 ) 

The integration over the rest of the variables is simplified, particularly in the context 
of the NLO corrections, by effecting a change of variables, namely 

(Y, r) — > « x° 2 ) = (v/7e y , ^ e~ Y ) (3.5) 
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Then it follows that 



^ J dz5(Y-Uog(^-^5(l-z)5(r-zx 1 x 2 ) \M ab ^'\ 2 T jf 



2tt 

= -q5$( x i -x\)5(x 2 -x%) 



\M ab ^ jj '\ 2 T jr 



(3.6) 



2 = 1 



rendering the remaining integrals trivial and thereby giving us the Born-level result 
for the ^-distribution. Having set the formalism, we may now calculate the next- 
to-leading-order contribution to the same. This involves the computation of matrix 
element squared for the processes a{k\) + &(A: 2 ) — > j(q) + c(k). The (1 + l)-body 
phase space integration can be performed in the CM frame of the incoming partons 
wherein the particle momenta may be parametrised as 

h = ^(i,o,. ..,o,i) 
h = _(i,o,-. -,o,-i) 

-q = —1 + 3,0,. •• 



-(l-z) sin0,-(l cos6) 

z, 0, • • • , (1 — z) sin 9, (1 — z) cos 6) 
Writing cos6> = 2y — 1, the two delta functions reduce to 



-»lr-ilog^ 
2 pi-q 



2 10g x 2 (z + y(l-z)) 



S(t — zxix 2 ) 



2xia;2a;ia;2(xia;2 + x\x^) 

{X\X2 — x\x2){xiX2 + x\x2) 2 
1 



X\ X2 



5(z-z*) 



where, 



x 2 x%{x 1 +x° 1 )(x 1 



(X1X2 — XiX2)(xiX 2 + X1X2) ' 

The above relations can be used to obtain 
J dPS 1+1 J dzS^Y-^ log 5(t - zx lX2 ) \M ab ^ 3j '\ 2 T jf 



™0™0 
x l x 2 

X1X2 



1 ( Q 



2\ e/2 



2x®X2(xiX 2 + 2^2) 



8tt V 4:n J T(l + e/2) x x x2(x x x\ + a+£ 2 ) 2 

• ( Xl _ x 0) (a . 2 _ 3.0) ^ + 3.0) (a . a + 3.0) X e/2 



\M ab ^ jj '\ 2 T jf 



(3.7) 



(3. 



J y=y*,z=z* 



(3.9) 



x 
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To obtain the contribution to the Y distribution from real gluon emissions, we sub- 
stitute eqn. (|3.9| ) in eqn.( p.2|) . Similarly, the virtual corrections can be obtained using 
eqn.( |3.4| ) with oneloop corrected matrix elements. The soft singularities cancel after 
adding the real emission contributions and virtual corrections to the Born process. 
The remaining collinear divergences are removed by mass factorization, or, in other 
words, by replacing the bare parton distribution with the renormalized ones using 
the Alteralli-Parisi kernels as follows 



/f(*) = £ r ^®#WS0 



(3.10) 



which implies 



ff(z) 



/ P 0,/4) 



€ 

e 

a s n f 



rS®L p (^/x|)+rg®/^ )M |) 



+r (i) 

' 99 



fg(z,^ F ) 



(3.11) 



Thus, we finally have, for the one-loop corrected distributions in the dilepton pair 
rapidity, 



£ F, SM+VA,q 
i=q 



D q q {^i, X®, fJ> F ) 



+ £ ? SP,q 
i=q 



(3.12) 



D„x (x®, x®, Hp) 



gq 



^qq{ X °i-i X 2) Vf) + D qq (x®, X®, Hp) 



with the constants FsM+VA,q an d J~sp, q as in eqn.( p.22"|) . The functions D can be 
split conveniently into the Born-approximation piece and the NLO corrections, viz 

DUx^xlti) = D^\xUlti) + a s D^\xlxlti) ( V = SM,VA,SP). 

(3.13) 

Once again, the analytical expressions for the \^A-type contact interactions are the 
same as those obtained within the SM and can be found in Ref. [26]. As for the 
S'P-type interactions, while the leading-order expression is simple 
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the NLO results are more complicated. Defining, for convenience, certain constants 



«oi = In 



K Cl = In 



2Q 2 (l-x° 2 ) { Xl -xl) 
H% (xi + x?) x% 
2x\ 
Xi + x\ 



n bl = In 



k 12 = In 



Q 2 (l-x° 2 )(x 1 -x° 1 ) 



I / 2 /y»0 ,yj 



(1 -*?)(!-:$ 



(3.15) 



o ^.o 



37 1 37 



1 ^2 



we have 



{x^i X2) l^ip) 



11 



2%C f 
N 



ii 



j dxiipf + J d Xl dx 2 (pf^ + (l ^ 2) 



H^xl, x° 2 , -2 + k( 2 + 6 C(2) + (3 + 2 k 12 ) In 



/4 



<pl q 



2K bl 



H qq; l(xi, X 2 , Hp) 



0u\ 37-^ 
+ Hgg(xi, X®, ftp) 



K 



ai 



+ 



2k c 



1 + K 







<p? = 



+ 



X\ 37 1 37-^ 

qq,u(Xl,X 2 , fJ, 2 p) X 2 + X 2 



ai x° 



(xi — Xi)(x 2 — x 2 ) (xi — Xi) X 2 



H q q( Xl ,X 2 ,fJ J 2 p) 



2,7^2 ™2 
12 



>i + a;?) (^2 + s°) + 



//qq,l(Xl,X 2 ,^) 

9 9 0^ 
12 ~~ ~ 1 2 



(xi + x?) (rr 2 + ^2) 

(3.16) 



and 



Dgq ^(^l) ^2' A*f) 



2ttT) 
TV 



X 



¥1 + J yp 2 - xi {X2+X ^ {X13 » +X2X ^ j 



(fif = H gq (x 1 ,xl,IJ, 2 p) ( 2 X \( X i - X \) + K ai [X° + (Xi - x\ 



„0\2 



91 



H g q j2 (x 1: X 2: H 2 p) 



xf + (xi — Xi 



0\2 



¥2 

x 2 - x 2 

— 37-^ *^"2p^'2 *^ 1 *^ 1 3/2*37 2 (3372 ~~ \~ ^^2) 



- \~ X -^X 3^23/2 ^^3^2 ~~ \~ 23^2 ^ ~~ \~ 237-^ 372 1*^2 ~~ 23/23/2 ~~ I - 23723/2 ~~ I - 23/2 ) 

H-237;[37^ 372 ( — ^2 ^2*^2 ^~ ^%2*^2 2^237^ ~~h 237^ ) 

1 2 3 / 5 _ 4 4 _ 4 3 2 , 9 2 3 , 9 4 , 9~0 5 \ 
I ju -^ju v 2 j: t.c 2 2 2 2 ~~ 2 2 ~~ ^ «x- 2 *■*-' 2 ~~ r~ 2 / 



with 



D^ix^xlul) = I% p >\4,4,ti)\(i~2) 



(3.17) 
(3.18) 
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and we have used the following notations, 



Hab,i 2 (xi, x 2 , = H ab (xi, x 2 , fi F ) - Hab(x1, x 2 , ftp) - F^xi, x°, fi F ) 
~\~H a t,(xi, x®, up) 

(3.19) 

H ab ,i{x 1 ,x 2 ,n F ) = H ab (xi,x 2 ,pr F ) - H ab (x° 1 ,x 2 ,fi F ) 

H a b,2(x 1 ,X 2 ,fI F ) = H ab (xi,X 2 ,[l F ) - H a b{xi,xl,H F ). 

4. Results and Discussion 

We now present numerical results relevant for the Run II of the Tevatron (y~S = 
1.96 TeV) as well as for the LHC (y/S = 14 TeV). Although our goal, namely the 
calculation of the NLO QCD corrections, would be quite independent of the value 
of the contact interaction scale A, for defmiteness we choose A = 6 (20) TeV for the 
Tevatron (LHC). Furthermore, in presenting our results, we shall consider only one 
of the couplings rf AB and ^ q AB to be non-zero and of unit strength. 
For the sake of convenience, we parametrize the cross section as 

o = °sm + 0"intf + o r] i (for the VA case) 

(4.1) 

a = cr SM + crg2 (for the SP case) 

and similarly for the differential cross sections. This has the advantage in that the 
total cross sections, for an arbitrary value of A can be easily reconstructed. 

4.1 The invariant mass distribution for the dilepton pair 

In Fig.|l|, we present the invariant mass distributions corresponding to the VA case. 
Note that the a v 2 piece (and, similarly, the a^. piece) depends only on the identity of 
the quark q taking part in the contact interaction and is independent of the chirality 
structure of the coupling. The interference term, on the other hand, does depend on 
the chirality structure as Fig.[l](a) amply demonstrates. As for <tsm, the rapid decrease 
in cross section with M is reflective of both the s _1 fall of the parton-level cross 
section as well as the rapid fall in parton distribution functions at higher momentum 
fractions. That the interference terms do not fall as fast is a consequence of the higher 
dimensional nature of the contact interaction Lagrangian. This, naturally, is even 
more evident for the 0^2 (cx^a) piece. Consequently, at high M values, the contact 
interaction contribution dominates over the SM piece. For the LHC, this dominance 
occurs at a larger M value as compared to the case of the Tevatron precisely because 
we have chosen to work with a much larger value of A for the former environment. 
And, expectedly, for identical couplings, the cross section due to a uu initial state 
dominates that originating from a dd initial state. In Fig.|l], we have chosen to limit 
ourselves only to these two initial states as the cross sections corresponding to the 
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Figure 1: The differential inclusive dilepton production cross-sections (at NLO) for the 
contact interaction terms. In each case, only one coupling (t],£) is assumed to be non- 
zero and of unit size. Also shown, for comparison, is the total SM contribution. The top 
and bottom panels refer to the Tevatron (\f~S~ = 1.96 TeV and A = 6 TeV) and the LHC 
(y/S = 14 TeV and A = 20 TeV) respectively. The right and left panels refer respectively 
to the pure contact interaction term and the interference with the SM. 



heavier quarks would be suppressed even further (note though that the experimental 
bounds on A is relaxed too for such cases). 

We should clarify, at this stage, that, in calculating the NLO cross sections 
shown in Fig.|I], we have made a specific choice of the renormalization scale and 
the factorization scale fip, namely, 

fi R = fi F = M= y/Q 1 . 

Postponing, until later, a discussion of the dependence on the scale choice, we may 
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Figure 2: The K -factors for the differential (in dilepton invariant mass) cross-section for 
(a) the Tevatron Run II and (b) the LHC. For the contact interactions, the K-factors are 
independent of the chirality structure of the operators, but depend on whether they are of 
the V A or the SP type. 

define now a invariant mass-dependent .fT-factor, namely 



'da LO (M)' 


-l 


~da NLO (My 


dM 




dM 



(4.2) 



where q refers to the identity of the quark, and the LO (NLO) cross sections are 
computed by convoluting the corresponding parton-level cross sections with the LO 
(NLO) parton distribution functions. In Fig.|2|, we exhibit the variation of Km with 
M for different choices of q. 

As derived in the previous section, and as already evinced in Fig.jl], the fractional 
correction depends only on the spin structure of the vertex, and not on the chirality. 
Thus, for a given quark, the i^-factor would depend on whether the interaction is of 
SP or V A type, but within each class, the chirality structure (namely whether it is 
LL, RR, LR or RL type) is quite irrelevant. The last statement also implies that, 
for the VA-type interaction, the i^-factor would be exactly the same as in the SM, 
as far as the particular quark initial state is concerned. Numerically, this feature is 
displayed in Fig.^. Of course, the i^-factor does depend on the identity of q. As can 
be expected, K M (VA) and K M (V A) are relatively close to each other and, in turn, 
to K M (SM). In fact, to a large measure, K M (SM) is but the weighted average of the 
other two, with the relative strengths being determined by the quark fluxes. That 
these .ff -factors fall monotonically with M for the case of the Tevatron and not so for 
the LHC is understandable in the light of the fact that, the former is a pp machine, 
while the latter is a pp one. As for K M (VA), the steep rise at large M values is 
but a reflection of the dominance of the Compton-like subprocess (sg — > £ + £~s and 
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sg — > f + frs) owing to the larger flux of gluons as compared to s/s, especially for 
large momentum fractions. 

The results for the S"P-type interactions are qualitatively similar, though quanti- 
tatively the .fT-factors are significantly larger than those for the VA-type interaction 
(or the SM). The numerical differences are but consequences of the the structures of 
the respective matrix elements. On closer inspection, K M (SP), for a given M, turns 
out to be the same as that for resonance production of a scalar /pseudoscalar of mass 
M [29]. 



4.2 The rapidity distributions 

We now turn to the distribution in a different kinematical variable, namely Y, the 
rapidity of the lepton pair 1 . However, rather than look at da/dY itself, we shall 
rather consider on d 2 a/dM dY, for this allows us to accentuate the effect of the 
contact interactions by concentrating on a suitable M range. 

At the LO, the variable Y is just a measure of the boost of the partonic center 
of mass with respect to the laboratory frame. At the NLO, one has to consider 
the effect of the initial state radiation as well. In either case, it is easy to see that 
d 2 a/dM dY (and, hence, da/dY) is an even function of Y. In Fig.[3], we exhibit the 
dependence of d 2 a/dM dY on Y, for a fixed value of M. As expected, for a large 
enough value of the latter, the effect of the contact interaction is clearly discernible 
and especially for the central rapidity region. Note that the contact interaction cross 
sections are significantly flatter in Y than the SM contribution. Once again, this is 
a reflection of the structure of the new physics matrix element as compared to that 
due to j/Z exchange. 

Analogous to K M defined in the previous subsection, one may now define a 
Independent X-factor of the form defined as 



Ky EE 



da LO (M,Y) 
dM dY 



-i -1 r 



da NLO (M,Y) 
dM dY 



(4.3) 



and we plot this quantity as a function of Y in Fig.f|. The results are quite reminiscent 
of those for Km (as displayed in Fig.|2]). It is noteworthy that, for the LHC, K Y 
shows a large upward swing at large Y, whereas K S M had seemed better behaved. 
The reason is not difficult to fathom. For the range of M spanned in Fig.|2](&), the 
cross section integral typically samples relatively moderate values of the Bjorken-x 
as compared to the higher- M« regime for the Tevatron case (Fig.Wa)). On the other 
hand, a phase space point such as (M = 700 GeV, |Y| = 2.5) necessarily pushes one 
to larger momentum fractions for the partons and thus, once again, it is the ratio of 
the strange-quark flux to that of the gluon that causes the upward turn in Fig.|(fc). 



lr This should be distinguished from the rapidity of an individual lepton. 
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Figure 3: As in FigQ but for the double differential (in rapidity Y and mass M) distri- 
bution instead. 



4.3 The choice of scale 



Until now, we have chosen each of the factorization scale \ip (relevant to both the 
LO as well as NLO calculations) and the renormalization scale fiR (relevant only for 
the NLO case) to be the same as the dilepton invariant mass M. As is well known, 
this choice is arbitrary and there is no theoretical guideline for making such a choice. 
Maintaining, for reasons of simplicity, /ir = fiF, we now examine the dependence of 
our calculations on this choice. To quantify the scale dependence of our result, we 
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Figure 4: As in Fig^, but for the double differential (in rapidity Y and mass M) distri- 
bution instead. 



define ratios Rm and Ry 
R Y M = 



dai(M, fx F = M) 
dM 

daj(M : Y, fj,p = M) 
dM dY 



dai(M, n F ) 
dM 

da I {M,Y J ii F ) 
dM dY 



(4.4) 



where I = LO, NLO. A value of R^y^f) close to unity would then signify a low 
sensitivity to the choice of scale and hence a more robust result. 

In Fig.|5], we display the above ratios for the case of the LHC and the V A inter- 
actions. Note that the variation of the cross section with the factorization scale is 
relatively small. Furthermore, the variation reduces significantly as one goes from 
the LO to the NLO case. This immediately points to the increased robustness of 
the prediction on inclusion of the corrections, and lends hope that the remaining 
scale ambiguity can, presumably, be reduced by adding still higher order corrections. 
Note that, at the leading order, these ratios are independent of the dynamics and 
reflect only the effect of the choice of the factorization scale on the parton densities. 
In other words, 



R L M U (SP) = R L A y(VA) } R Y U (SP) = R Y U (VA) . 

At the next-to-leading order, the dynamics does play a role. However, the differences 
between the i?-ratios for the SP and V A cases are too small to be noticeable on the 
scale of Fig.[5[ The results are similar for the case of the Tevatron as well. 
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Figure 5: The K -factors for the differential (in dilepton invariant mass) cross-section for 
(a) the Tevatron Run II and (b) the LHC. For the contact interactions, the K-factors are 
independent of the chirality structure of the operators, but depend on whether they are of 
the V A or the SP type. 

5. Conclusions 

To summarize, we have performed a systematic calculation of the next-to-leading 
order QCD corrections for the Drell-Yan process in theories with contact interactions. 
Contrary to naive expectations, we demonstrate explicitly that the QCD corrections 
are meaningful and reliable in the sense that no undetermined parameters need be 
introduced. 

We have analyzed both the invariant mass distribution and the rapidity distribu- 
tions for the dilepton pair at either of the Tevatron and the LHC. The enhancements 
over the LO expectations are found to be quite significant. The corresponding K- 
factors are presented in a form suitable for use in experimental analyzes. 

For the VA-type interactions, the analytical structure of the corrections are sim- 
ilar to those for the SM. However, a significant dependence on the flavour structure is 
found and needs to be carefully accounted for in obtaining any experimental bounds. 
For the SP-type interaction, not only are the analytical results quite different, but 
the consequent .fT-factors are typically larger than those within the SM. 

Finally, we have investigated the sensitivity of our results to both the factoriza- 
tion and renormalization scales. As expected, we find such dependences to be greatly 
reduced for the case of the NLO results as compared to that for the LO case. This 
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indicates the robustness of the calculations. 
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